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Abstract 



We address the issue of generating cutting planes for mixed integer programs from 
multiple rows of the simplex tableau with the tools of disjunctive programming. A 
cut from q rows of the simplex tableau is an intersection cuts from a g-dimensional 
parametric cross-polytope, which can also be viewed as a disjunctive cut from a 2 9 -term 
disjunction. We define the disjunctive hull of the q-row problem, describe its relation 
to the integer hull, and show how to generate its facets. For the case of binary basic 
variables, we derive cuts from the stronger disjunctions whose terms are equations. 
We give cut strengthening procedures using the integrality of the nonbasic variables 
for both the integer and the binary case. Finally, we discuss some computational 
CO ' experiments. 

\o ' 

1 Introduction: intersection cuts and disjunc- 



tive programming 



In the last few years a considerable effort has been devoted to generating valid cuts for 
^ \ mixed integer programs from multiple rows of the simplex tableau, with a focus on cuts 

from two rows. This research was pioneered by the 2007 paper of Andersen, Louveaux, 
Weismantel and Wolsey PQ, followed by Borozan and Cornuejols [13J, Cornuejols and 
Margot [T7], Dey and Wolsey [19] and many others (|1H [12, 18, 20] ; for a recent survey 
see P]). 

All of these papers view and derive the multiple-row cuts as intersection cuts, a 
concept introduced in [2], i.e. cuts obtained by intersecting the extreme rays of the 
cone defined by a basic linear programming solution with the boundary of a convex 
set whose interior contains no feasible integer point. Intersection cuts are equivalent 
to disjunctive cuts, and in this paper we apply the tools of disjunctive programming 
to the study of cuts from multiple rows of the simplex tableau. Two early versions of 
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this paper were presented at the 2009 Spring Meeting of the AMS in San Francisco [5] 
and at the 20 th ISMP in Chicago [6]. 

The structure of our paper is as follows. In the remainder of this section we out- 
line the connection of intersection cuts with disjunctive programming. In section [2] we 
introduce the concept of disjunctive hull associated with q rows of the simplex tableau 
and examine the relation between the disjunctive hull and the integer hull. We then 
give a geometric interpretation of cuts from q rows of the simplex tableau as cuts from 
a g-dimensional parametric cross-polytope (section [3]) , followed by a theorem relating 
the facets of the disjunctive hull to those of the integer hull (section [3]) . In section [5] we 
specialize these results to the case of q = 2. The next section (JH]) discusses the strength- 
ening of our cuts when some of the nonbasic variables are integer-constrained. Section [7] 
deals with the 0-1 case, when the stronger disjunction whose terms are equations can 
be used to derive stronger cuts. Finally, section [8] describes some computational ex- 
periments. 



* 



Suppose a Mixed Integer Program is given in the form of q rows of the simplex 
tableau 

x = x + ^2r j Sj , xeZ^sERl (1.1) 
ieJ 

where x is a basic feasible solution to LP, the linear programming relaxation of a MIP, 
and we are interested in generating an inequality that cuts off x but no feasible integer 
point. 

Theorem 1.1. (Balas ]2$)- Let TCI' be a closed convex set whose interior contains 
x but no feasible integer point. For j G J, let s* := max{sj : x + r 3 Sj G T}. Then the 
inequality as > 1, where ay = 4f ; j G J, cuts off x but no feasible integer point. 

The inequality as > 1 is known as an intersection cut. Theorem 11.11 is illustrated 




by Figure [TJ In both cases (a) and (b) the convex set T consists of the intersection 
of two halfspaces, but in (b) the two halfspaces are defined by hyperplanes parallel 
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to one of the coordinate axes, and so their intersection defines an infinite strip. The 
intersection cut from this latter set T is the Gomory Mixed Integer cut (GMI) |21j . 

This particular class of intersection cuts, the GMI cuts, has played a crucial role 
in making mixed integer programs practically solvable. These cuts are derived from a 
convex set of the form [xi\ < Xi < \xi~\, where Xi = xi + Y^jej r } s j * s one °f the rows 
of an optimal simplex tableau and [xi\ < x~i < \xi\. More generally, cuts obtained 
from a convex set of the form itq < irx < itq + 1, where (tt, ttq) is an integer vector 
with gcd(7r) = 1, are known in the literature as split cuts [16] , It is then natural to ask 
the question whether intersection cuts derived simultaneously from several rows of a 
simplex tableau have some properties that distinguish them from split cuts. It was this 
question that has led to the investigation of intersection cuts from maximal lattice-free 
convex sets by [TJ [13] and others. 

We propose a different approach to the same problem, which promises some com- 
putational advantages. The approach is that of Disjunctive Programming, a nat- 
ural outgrowth of the study of intersection cuts. To see the connection, consider 
an intersection cut from a polyhedral set with the required properties, of the form 
T := {x : d l x < d l , i = 1, . . . ,m}. Clearly, the requirement that intT should con- 
tain no feasible integer point, can be rephrased as the requirement that every feasible 
integer point should satisfy at least one of the weak complements of the inequalities 
defining T, i.e. should satisfy the disjunction 

m 

V(^>4)- (i-2) 
i=i 

Therefore an intersection cut from T can be viewed as a disjunctive cut from (jl.2p . 
While these two cuts are essentially the same, the disjunctive point of view opens up 
new perspectives. Thus, suppose that in addition to (|1.2p . all feasible solutions have 
to satisfy the inequalities Ax > b. Then one way to proceed is to generate all valid 
cutting planes from (jl . 2[) and append these to Ax > b. The resulting system will be 

P := i^x € M. n : (Ax > b) H conv 

But another way to proceed is to introduce Ax > b into each term of the disjunction 
(fL2j) . i.e. replace dH2J) with 

V ( tJ_\ ) ■ ™ 

and take the convex hull of this union of polyhedra: 

Q:= ^(£(^14)) 

Now it is not hard to see that Q C P, and in fact Q is in most cases a much tighter 
constraint set than P. We illustrate the difference on a 2-term disjunction. Given an 
arbitrary Mixed Integer Program, let (tt, ttq) be an integer vector with a component TTj 
for every integer-constrained variable. Then the disjunctive cut derived from 

TTX < TTq V 1TX > 7TQ + 1 (1-4) 
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is equivalent to the intersection cut derived from the convex set 

7T() < TTX < TTq + 1, 

illustrated in Figure [TJ On the other hand, the disjunction 

I AX ~ b \ V ( AX ~ b \ (15) 

\ TTX < TTq J \ TTX > TTq + 1 J 

gives rise to an entire family of cuts, whose members are determined by the multipliers 
u, v associated with Ax > b in the two terms of this more general disjunction 

(-7T — uA)x < TTq — ub V (jT + vA)x > TTq + vb + 1 (1-6) 

Cuts derived from a disjunction of the form (|1.4p are called split cuts, a term that 
reflects the fact that (11.40 splits the space into two disjoint half-spaces. Cook, Kannan 
and Schrijver |16j who coined this term also extended it to the much larger family of 
cuts derived from disjunctions of the form (|1.6p . 

Disjunctive sets of the form (jl.3p or (jl.5p represent unions of polyhedra, and the 
study of optimization over unions of polyhedra is known as Disjunctive Programming. 
Its two basic results are a compact representation of the convex hull of a union of poly- 
hedra in a higher dimensional space, and the sequential convexifiability of facial dis- 
junctive sets [HE]. The application of disjunctive programming to mixed 0-1 programs 
has become known as the lift-and-project method [7]. Here we apply this approach to 
the study of intersection cuts from multiple rows of the simplex tableau. 



2 Integer and disjunctive hulls 

Consider again a system defined by q rows of the simplex tableau, this time without 
the integrality constraints: 

x = f + J2r j sj, sj>0, j £ J, (2.1) 

where /, r J £ W, j £ J := {1, . . . ,n}, and assume < fa < 1, i £ Q := {1, . . . , q}. 
This assumption can be made without loss of generality since setting x\ = xi — [fi\ 
and f- = fi — L/iJ , i £ Q, we have that x[, f-, i £ Q satisfy the assumption. The set 

P L ■= s ) G M.1 x R n : (x, s) satisfies flUE])} (2.2) 

is the polyhedral cone with apex at (x, s) = (/, 0) defined by the constraints that are 
tight for this particular basic solution. Imposing the integrality constraints on the basic 
components we get the mixed integer set 

Pj := {(x,s) £ Pl : Xi integer, i £ Q}, (2.3) 

whose convex hull, convP/, is Gomory's corner polyhedron [22j, or the integer hull 
of the MIP over the cone Pl. The main objective of the papers mentioned in the 
introduction was to study the structure of Pj for small q, with a view of characterizing 
the facets of conv P/ and minimal valid inequalities for Pj . 
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Consider now the following disjunctive relaxation of Pj, obtained by replacing the 
integrality constraints on X{ with the simple disjunctions X{ < V X{ > 1, i G Q: 



P D := {(x, s) G P L : Xi < V a* > 1, iGQ}. 



(2.4) 



Like Pj, Pd is a nonconvex set. Its convex hull, convPo, which we call the simple 
disjunctive hull, is a weaker relaxation of Pj than convP/, i.e. convPo D conv Pj, but 
it is easier to handle, since it is the convex hull of the union of 2 q polyhedra. Thus one 
can apply disjunctive programming and lift-and-project techniques to generate facets 
of convPo at a computational cost that for small q seems acceptable. In this context, 
the crucial question is of course, how much weaker is the relaxation convPo than 
convPf? We will pose this question in a more specific form that will enable us to give 
it a practically useful answer: when is it that a facet defining inequality for convPo is 
also facet defining for convP/? In other words, which facets of the (simple) disjunctive 
hull are also facets of the integer hull? Before addressing this question, however, we will 
take a side-step, by introducing a third kind of hull. If we strengthen the disjunctive 
relaxation of Pj by replacing the inequalities in the disjunctions x^ < V X{ > 1, 
i £ Q, with equations, we get the set 



Pp := {(x, s) G P L : x t = V x t = 1, i G Q}, 



(2.5) 



whose convex hull, convP^, we call the 0-1 disjunctive hull. For a general mixed 
integer program, the 0-1 Disjunctive Hull is not a valid relaxation, in that it may cut 
off nonbinary feasible integer points. Indeed, we have 



conv 



Pd 12 conv Pj D conv P D , 



where both inclusions are strict and are valid in the context of mixed integer 0-1 
programs only, since all the non-0-1 integer points that it cuts off are infeasible. Hence 
convP^ is equivalent to the convex hull of Pj PI {x : Xi < 1, i£ Q}, or the integer hull 
of Pj reinforced with the bounds on the X{. However, as we will see later on, finding 
facets of convP^ requires roughly the same computational effort as finding facets of 
conv Pjj. 

The upshot of this is that for the important class of mixed integer 0-1 programs, 
all facet defining inequalities of convP^ are facet defining for the integer hull. Fur- 
thermore, from the sequential convexification theorem of disjunctive programming, all 
such inequalities are of split rank < q, i.e. they can be obtained by applying a split 
cut generating procedure at most q times recursively. 

The set Pd of f)2.4[) is the collection of those points (x, s) G M. q x R n satisfying f)2.1|) 
and Xi < V Xi > 1, i G Q. Put in disjunctive normal form, this last constraint set 
becomes 

/ xt<0\ ( x x > 1 \ / xi > 1 \ 



x 2 < 



V x q < J 



V 



x 2 < 



V x q < J 



V 



V 



x 2 > 1 



V x q > 1 J 



(2.6) 



Each term of (12. 6p defines an orthant-cone with apex at a vertex of the (/-dimensional 
unit cube. These 2 q orthant-cones are illustrated for q = 2 in Figure [2j 
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1,1 




Using (|2.ip to eliminate the ^-components and denoting by r l the i-th. row of the 
q x n matrix R = (r- ? )" =1 , (j2.6j) can be represented in R n as s > and 



-r s 



> /i\ 



-r 2 s > / 2 



V 



/ rh > 1 - /i \ 

-r 2 s > / 2 



V 



V 



/ 



( r x s > 1 - h \ 
r 2 s > l-/ 2 

V r«a > 1 - /, J 



(2.7) 



If P (n) C 



denotes the polyhedron defined by the i-th term of this disjunction 



plus the constraints s > 0, then can be defined in n-space as Pj 
where t = 2 q . Furthermore, we have the following: 



(■») 



D 



u 



i=l J 



(n) 

Theorem 2.1. convP^, is the set of those s 6 R n satisfying s > anc? aZZ f/ie 
inequalities as > 1 whose coefficient vectors a € R n satisfy the system 



a + r 1 tiii + ■ ■ ■ + r 9 n 



la 



- r 1 U2i + • ■ ■ + r 9 M2ij 



> 

> 



/i«i 



' + fqU\q 



(l-/l)"21 H 1" / 9 M2q 



- r 1 «ii — ■ ■ ■ — r q ut q > 

> 1 

> 1 



(2.8) 



for some Uik > 0, i = 1, . . . ,t = 2 q , k = 1, . . . , q. 



{i-h)un + ■■■ + (!- f q )u tq >l 



(n) 

Proof. Applying the basic theorem of Disjunctive Programming to convP^, we in- 
troduce auxiliary variables s 1 £ R n , 6 R, i = 1, . . . , i = 2 9 , and obtain the higher- 
dimensional representation 
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-r 2 s 1 



-fizi 
-hzi 



= 

> o 

> 



r l g 2 



-r 2 s 2 



-fqZl 



> o 

> o 

> 



- r 1 s 2 



~fqZ2 



> 



(2.9) 



2 f 

-r s 



-(l-fi)z t > 

-(1-/2)* > 



-Q--fq)Zt > 

z\ +z 2 H = 1 



s l > 0, i = Zi > 0, i = 

Projecting this system onto the s-space with multipliers a; uu, . . . , U21, • • • , U2 q ; 
; Uti , . . . , Ut q , we obtain 



o + r 1 ^!! +■■■+ 



> 



-riut q > 



-/3 + /iun H h 



> (2.10) 



-/8 +(l-/i)«t/i +■•■+ (1-/,S > 

u ik > 0, i = 1, • ■ • ,i, k = 1, . . . , q 

Applying the normalization (3 = 1 (clearly /3 = — 1 does not yield any cuts since it 
makes (|2.10p unbounded) we obtain the representation given in the theorem. □ 

In order to restate the system (|2,8p in a more concise form, for each i € {1, . . . , t} 
we partition the index set Q := {1, . . . , q} into 

Qf := {k G Q : Uik has coefficient vector r^} 
:= {k £ Q : Ujfc has coefficient vector — r^}, 
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with Qj U QT = Q, i = 1, . . . , t = 2 k . Then (|2l!j) can be restated as 

a + Yj ( r ^ :k£ Qt) Yl { rkuik :k ^Qi) ^ 

YifkUik ■ k € Q+) + £((1 - fk)uik :keQ7) > 1, i = l,...,t D) 

^jfc > 0, i = 1, . . . , t = 2 q , k = l,...,q 

The system f|2.8t ) has several interesting properties described in the next few propo- 
sitions. 

Proposition 2.2. For any p € W 1 , p > 0, all optimal basic solutions to the cut 
generating linear program 

mm{pa : (a,u) satisfies f2.8\ )} (CGLP)q 

are of the form 

ctj = maxjaj, . . . , a*}, (2-H) 

where 

a i : = - I>* u ^ '■ k G Qt) + '■ k G QD» ( 2 - 12 ) 

i = 1, . . . ,t = 2 q , with the Uik satisfying t2.8\ ' ). 
Proof. The constraints of (|2.8I ) require 

ay > « = 1, • • • ,t, j = 1,. • • ,n 

Suppose there is an optimal solution to (CGLP)q such that ay, > maxja^ : i = 
1, ... ,i} for some j* 6 {1, . . . , n}. Then setting ay t equal to the maximum on the 
righthand side, and leaving ctj unchanged for all j ^ j* yields a better solution, contrary 
to the assumption. □ 

in) 

Proposition 2.3. In any valid inequality as > 1 for convi-^) > a j > 0, j = 1, . . . , n. 

Proof. From (|2.1ip . aj > a*- for all i = 1, . . . , 2 9 , and in view of the presence of all 
sign patterns of r^Uik in the expressions (j2. 12|) . there is always an index i £ {i, . . . , 2 9 } 
with a) > 0. □ 

Proposition 2.4. For any frasic solution (a,u) to (CGLP)q that satisfies as strict 
inequality some of the nonhomogeneous constraints of A2.8\ ' ). there exists a basic solu- 
tion (a,u), with a = a, that satisfies at equality all the nonhomogeneous constraints of 
(CGLP) Q . 

Proof. Let (a, u) be a basic solution to (CGLP)q that satisfies as strict inequality some 
of the nonhomogeneous constraints of (12.81 ). W.l.o.g., assume that 

fiun H h f q ui q -6 = 1 

is one of those constraints with the surplus variable 6 positive in the solution (a, u). 
We will show that there exists a solution (a, u), with a = a and Uik = Uik for alii ^ 1 
and all k, such that 

flUn H h f q Ui q = 1. 
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Applying this argument recursively then proves the Proposition. 

Fix all variables of (CGLP)q except for tin, ■ • • jU\ q , at their values in the current 
solution. The fixing includes all the surplus variables except those in the n + 1 rows 
containing un, . . . ,ui q . This leaves the following constraint set in the free variables: 

-rjun r q ju lq +tj = aj j = 1, . . . , n 

+ — + f g ui q -9 = 1 (2.13) 
un, ■■-,u lg >0, tj > 0, j = 1, ... ,n, 9 > 

Here 6, tj represent the surplus variables of the respective constraints. We claim that 
this system has a solution with 9 = 0. To see this, consider the linear program 

min{6> : Uik,tj and 9 satisfy (|2.13p | 

and its dual, 

n 

max Ao + aj Aj 
3=1 

subject to 

/1A0 - E r}\ 3 < 

< 

< 1 

Aj < 0, j = l,...,n 

Since aj > 0, j = 1, . . . , n, it is not hard to see that the dual linear program has an 
optimal solution Ao = 0, Aj = 0, j = 1, ... ,n and hence the primal has an optimal 
solution with 9 = 0. □ 

The obvious and important consequence of Proposition 12.41 is that for all practical 
purposes we can replace all 2 q nonhomogeneous inequalities in the constraint set ([2. 81 ) 
of (CGLP)q with equations. In view of Proposition 12.21 it then follows that we may 
restrict our attention to basic feasible solutions that satisfy at equality n + 2 q out of 
the n x 2 q + 2 q inequalities of (|2.8I ) other than the nonnegativity constraints. 

At this point we introduce the characterization of convP^, the 0-1 disjunctive hull 
defined by (12. 5j) , closely related to that of conv Pp . Just as in the case of Pjj , we denote 
by P= (n) the union of polyhedra in R n representing the disjunction (|2.7[) in which all 
the inequalities have been replaced by equations. The following Theorem is the analog 
of Theorem 12.11 for this case. 

=(n) 

Theorem 2.5. convP^ is the set of those s € M. n satisfying s > and all inequalities 



n 

/gA - E r j Aj 
-A 
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as > f3 whose coefficients satisfy the system 



a + r l U\\ H h r l ui q 



> 







a 



— . . . — r <? 



wt 9 



> 








(2.14) 



P + fiuu H h 



-/3 



+ (l-/l)«fl + ' ' ' + iX~fq)utq 







/or some u^, i = 1, . . . , i = 2 9 , k = 1, . . . ,q. 

Proof. The proof of Theorem 12.11 goes through with the following modifications. Since 
the inequalities in the disjunction (j2.7|) are all replaced with equations, the inequalities 
in the system (|2.9|) . other than the nonnegativity constraints, also become equations. 
As a consequence, the variables Uik of the projected system (|2.10p become unrestricted 
in sign. The remaining difference between (|2.14|) and (|2.8p is the fact that in (|2.14|) the 
last 2 q constraints are equations rather than inequalities. This is due to the fact that 
Proposition 12.41 applies here too. In other words, if we denote by (|2.14l / ) the system 
obtained from (|2.14p by replacing the equations containing f3 with inequalities >, then 
for any basic solution (a, u) to (CGLP)q that satisfies as strict inequalities some of 
the constraints ()2.14l / ) containing /3, there exists a basic solution (a,u), with a = a, 
that satisfies at equality all the constraints containing j3. The proof is essentially the 
same as that of Proposition 12.41 

Thus the two basic differences between the systems describing conv P^ and conv P^ 
are that (a) the latter also contains inequalities of the form ax < 1 (corresponding to 
/3 < 0), and (b) the coefficients aj of the latter can be of any sign. □ 

We now return to the simple disjunctive hull, convPo, and describe its vertices. 

Proposition 2.6. Every vertex of conv P^ is a vertex of some P^ n \ i € {1, . . . ,2 9 }. 

Proof. Let v be a vertex of convPrj . If v E P^ for some i G {l,...,t = 2 q }, then 

(ri) 

v must be a vertex of P- , or else it could be expressed as a convex combination of 
points in , hence of P^ . On the other hand, if v ^ UP^ but v € convP^ , then v 
is a convex combination of points in \j\ =l P^ n \ hence of conv P^\ a contradiction. □ 

Next we describe the vertices of Pj , i € {1,...,2 9 }. We will call a vertex of 
conv P^ (of Pj ) integer if it defines an integer x through f)2 . 1 1) : in other words if 
fi + r l s is integer for i = 1, . . . , q. All other vertices will be called fractional. 

For any particular i* € {1, • • • , 2 9 }, 



where (Qi t , Q \ Qi t ) is the partition of Q that defines i*. 

(n) 

Proposition 2.7. P/ L can have three kinds of vertices, distinguished by the corre- 
sponding x-vectors that belong to one of these types: 

(a) 0-1 vertices: Xh = 0, h e Qi t and Xh = 1, h G Q \ Qi f . 



P£> := {s € : r h s < -f h , h G Q u , r h s >l - f h , h G Q\ Q t J 
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(b) non-binary integer vertices: Xh € Z_, h € Qi t , Xh G Z+, h £ Q \ (here Z_ 
and Z + stand for the nonpositive and nonnegative integers respectively). 

(c) fractional vertices: n < 0, /j £ Qu> %h > 1> h E Q \ Qi. t , with at least one 
inequality strict. 

Proof. The three cases become exhaustive if the following fourth one is added: (d) frac- 
tional vertices with < Xh < 1 for some h € Q. But this case clearly violates at least 

(n) 

one of the constraints defining P- . □ 

(n) 

Note that -Fy can have several distinct vertices with the same associated x-vector, 
corresponding to basic solutions with the same x-component. Note also that if a 
component x^ of a vertex is fractional, then x^ < or x^ > 1. 

The next theorem characterizes the facets of the simple disjunctive hull. 

Theorem 2.8. The inequality as > 1 defines a facet of conv if and only if there 
exists an objective function of the linear program ( CGLP)q of Proposition \2.$\ with 
p > such that all optimal solutions (a,u) have a = a. 

Proof outline. This is a special case of Theorem 4.6 of [3j. The inequality ax > 1 
defines a facet of conv P^ if and only if a is a vertex of the polar of conv P^ , which 
is the projection of (|2.8p onto the a-space. But a is a vertex of this polar if and only if 
there exists an objective function vector p > such that pa attains its unique minimum 
at a. □ 

If the system (|2.4p defining Pl is of full row rank q, then the dimension of conv Pq 
is n, since there are q + n variables and q independent equations. The dimension of 
conv P^ is also n, so the facets of convPp are of dimension n — 1. 

From a computational standpoint, the most important feature of (CGLP)q is that 

(n) 

the facets of the n-dimensional convP^, can be generated by solving a smaller CGLP 
in a subspace of at most t = 2 q variables s,- , and lifting the resulting inequality into the 
full space. The idea of generating cuts in a subspace of the original higher dimensional 
cut generating linear program and then lifting them to the full space goes back to [10] , 
where lift-and-project cuts were generated from a 2-term disjunction by working in the 
subspace of the fractional variables of the LP solution. Here we are working with a 2 q - 
term disjunction, and are considering a different subspace, suggested by the structure 
of the problem at hand, but the lifting procedure is essentially the same as the one 
used in [3(8]. 

Since our cuts are derived from a disjunction with 2 q terms, if we want to create a 
subproblem in which all terms are represented, we need 2 q out of the n variables ay of 
our (CGLP)q. Furthermore, the 2 q vectors r 3 corresponding to these ay have to span 
the subspace M q of the x-variables. Solving the (CGLP)q in this subspace yields 2 q 
values ay and q x 2 q associated multipliers Uik, i = 1, • • • , 2 q , k = 1, . . . , q; and these 
multipliers can then be used to compute the remaining components of a, given by the 
expressions (I2.1ip and (|2.12j) . The significance of this is that the computational cost 
of generating facets of convPo grows only linearly with n. Of course this cost grows 
exponentially with q, but the approach discussed here is being considered for small q. 

The choice of the subspace is intimately related to the question of deciding which 
facets of the disjunctive hull are also facets of the integer hull. The best way to address 
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this question and that of the subspace to be chosen, is to first interpret the inequalities 
defining the disjunctive hull as intersection cuts. 



3 Geometric interpretation: Cuts from the q- 
dimensional parametric cross-polytope 

Consider the q-dimensional unit cube centered at (0, ... ,0), K q := {x € W 1 : — | < 
Xj < |, J € Q}- Its polar, K° := {x € M. q : xy < 1, Vx € K}, is known to be the 
^-dimensional octahedron or cross-polytope; which, when scaled so as to circumscribe 
the unit cube, is the outer polar of K q : 

K* = {x £ R q : \x\ < \q}, 

where \x\ = XXI x jl : j = 1, . . . ,q}. Equivalently, \x\ < \q can be written as the system 

— x\ — ■ ■ ■ — x„ < 



X\ — • • • — X q < 



Xl + • • • + Xq < 



k 
h 

(3.1) 

1 „ 



of t = 2 q inequalities in q variables. 

Moving the center of the coordinate system to , \) changes the righthand 

side coefficient of the i-th inequality in (13. ID from to a value equal to the sum of 
positive coefficients on the lefthand side of the inequality. Indeed, if q + and q~ denotes 
the number of positive and negative coefficients, then ^q + \q + — \q~ = q + . 

Next we introduce the parameters v ik, i = 1, • • • , t = 2 q , k = 1, . . . , g, to obtain the 
system 

-V U Xl - ■■■ - Vl q X q < 



V 2 lXl - ■■■ - V 2 qXg < V21 
-V31X1 + ••• - V 3q Xq < V 31 

VtlXl + • • • + VtqXq < V t \ + . . . + V tq 

Vik>0, i = l,...,t = 2 q , k = l,...,q. 
Note that the constraints of (13.21) are of the form 



(3.2) 



vikXk - ^2 vikXk - Vik -- 

k£Q+ keQ~ keQ+ 

where Qf and are the sets of indices for which the coefficient of x& is +vn~ and — Vik, 
respectively. Note also that all inequalities that have the same number of coefficients 
with the plus sign have the same righthand side, equal to the sum of these coefficients. 

The system (|3.2p is homogeneous in the parameters v^k, so every one of its inequal- 
ities can be normalized. Since we are looking for a connection with the system (I2.8P 
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defining (CGLP)q, we will use the normalization given by this system and Proposi- 
tion E31 i.e. 

/if 11 + • ■ ■ + fqVlg = 1 

(l-/l)«21 + •■■ + fqV 2 q = 1 

(1-/1M1 + ••• + (i-/,K = i 

Note that these normalization constraints are of the general form 

(i - /fc)wtA + e f kVik = l - 

h&Q+ h£Q- 

Let ET*(v) denote the parametric cross-polytope defined by (|3.2p and (|3.3p . It is not 
hard to see that for any fixed set of vn~, f)3.2|) defines a convex polyhedron in x-space 
that contains in its boundary all i£R' such that Xk € {0, 1}, k € Q, hence is suitable 
for generating intersection cuts. Furthermore, letting K*^ n \v) be the expression for 
K*(v) in the space of the s-variables, obtained by substituting / + Rs for x into (|3.2p . 
we have 



Theorem 3.1. For any values of the parameters Vik satisfying \3. £)) and \3. 3\) . the 

intersection cut as > 1 from K*( n \v) has coefficients a>j = where 

s* = max{sj : f + r j Sj E K< n \v)}. (3.4) 
Proof. This is a special case of Theorem 11.11 □ 

In order to compare the intersection cut as > 1 with the cut as > 1 from the 
(/-term disjunction (|2.7p . we have to restate (|3.4p in terms of the system of inequalities 
defining K< n \v). This means that / + r 3 s*j has to be expressed as the intersection 
point of the ray / + r J Sj, Sj > 0, with the first facet of K*( n \v) encountered. This 
yields 

s* = min{s], . . . , s*}, (3.5) 

where the s*- are obtained by substituting + Ylh=i r j s h f° r x k-> k = 1, . . . ,q into the 
z-th inequality of (j3.2j) . and setting s^ = for all h ^ j: 



s l j = max < 



( \ 

E Uifcr i ~ E VikT ) \ s i- E v ' ik ^ ~ fk)+ E VikJ * k 
\keQ+ keQj- J keQ+ k&Q- 



i = l,...,t = 2 q . 

Clearly, this maximum is bounded whenever the coefficient of Sj is positive, in which 
case, if we normalize by setting YlkeQ + Vik ^ ~ / fe ) + ^keQ~ v ' ik ^ k = 1 ' we obtain 

4 = ( E «*r? - E «*r? I ■ ( 3 - 6 ) 



13 



Comparing f)3.5[) and f|3.6[) to the expressions (|2.1ip and (|2.12p for the coefficient 
a,j of the lift-and-project cut as > 1 of Proposition 12.21 we find that setting = 
for all i, k, as well as Qf = and = Qf , we obtain 6y = ay. 

This proves 

Corollary 3.2. The intersection cut as > 1 from the parametric octahedron K*( n '(v) 
is the same as the lift-and-project cut as > 1 corresponding to the ( CGLP)q solution 
(a, u), with v ik = u ik , i = 1,. . . ,t, k = 1, . . . , q. 

4 Facets of the disjunctive hull and the integer 
hull 

Consider again the disjunctive relaxation of Pi 

P D = s ) G R q x R n : x = f + Rs, s > 0, x { < V x { > 1, i G Q} 

introduced at the beginning of section [21 where x,f G R q , R G R 9Xn , and Q := 
{1, . . . , q}. For i = 1, . . . , t = 2 q , let p % be the vertex of K q , the q-dimensional unit 
cube, defined by p\ = 0, i G Q 4 + , p k = 1, i £ Q^. 

Next we give a sufficient condition for an inequality as > 1 valid for to define 
a facet of convP/, which for small q leads to an efficient procedure for generating 
inequalities that are facet defining for conv Pj . 

(n) (n) 

The dimension of Pj being n > 2 q , as > 1 defines a facet of conv Pf' if there 
exists a subspace M 29 of W 1 such that the restriction of as > 1 to this subspace defines 
a facet of conv Pj ' . If this is the case, then the inequality in question can be lifted to 

(n) 

the full space to yield a facet of conv Pj by using the it-components of the solution 
(a, u) to the CGLP in the subspace to compute the missing coefficients ay. 

Theorem 4.1. Let as > 1 be a valid inequality for Pd corresponding to a basic solution 
(a,n) of (CGLP)q, and let p % , i = 1, . . . ,2 q , be the vertices of K q . Suppose for each 
p l , i = l,...,2 q , there exists a subset Jj C J containing the indices of q linearly 
independent rays r J1 , . . . , r Jq , and a vector A G M^_, satisfying 

P i -f=f2i^X j , f> = L (4.1) 

j=ji 3=3i 

Then the inequality Ylje.j a j s j — ^ defines a facet o/convPj' J '\ and its lifting based 
on the u-components of the solution (a,u) defines a facet of conv Pr . 

Proof. Suppose the subset of 2 q rays indexed by J satisfies the requirements of the 
Theorem. Then for every % = 1, . . . , 2 q , the vertex p % of K q satisfies 

p* = E(/-^')^ £a, = i 

3=31 3=31 

for some Xj > 0, j = j\, . . . ,j q , i.e. p % can be expressed as a convex combination of the 
q points /+— r J , j = j\, . . . , j g . But / + — r 3 = f + r^s* is the intersection point of the 
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ray / + r J Sj with bdK* hence each of these points satisfies as = 1 and consequently 

so does p i . Since Yl a j s j > 1 is satisfied at equality by 2 q integer points of convPj' J '\ 
jeJ 

it defines a facet of the latter. Furthermore, lifting the remaining coefficients aj of 
the inequality by using the u-components of (a, u) yields a facet defining inequality for 
conv/^ . U 

The sufficient condition of Theorem 14.11 is not necessary. There are two kinds of 
situations not satisfying the above condition, in which a valid inequality as > 1 for Pjj 
may define a facet of convP/. The first one involves an inequality as > 1 such that 
although (14. ip is not satisfied for all 2 q vertices p l of K g , nevertheless convPo has 2 q 
vertices whose x-components p l satisfy (|4.ip . i.e. convPo has multiple vertices with 
the same a>component. The second situation involves facet defining split cuts. 



5 The two-row case 



We now restrict our attention to the case q = 2, i.e. we consider two rows from a 
simplex tableau of a MIP problem with the variables xx,X2 and Sj,j E J: 



P L = {(x,s) €M 2+ I J I : X! =/i + E 

X2 =/ 2 + E 

Sj>o jeJ}. 



(5.1) 



where x\,X2 are basic variables required to be integers and Sj,j G J are non-basic. 
This is the case studied by Anderson, Louveaux, Weismantel and Wolsey [I]. Let 
Pi = {(x,s) € I? x Rl J l : (x, s) £ Pl}, and < /i,/2 < 1. The column vectors 
rj,j € J, represent the extreme rays of the cone in Rl^l with apex at (f\, /2). 

We will say that a ray r,- in ()5.ip hits an orthant-cone Qi,i € {1, ... ,4} if there 
exists Ao > such that / + Xrj G Q j for all A > Ao . 

For the case of 2 rows the disjunction (I2.7P becomes 



-r 1 * > /i' 
-r 2 ,>/ 2 .' 



r x s > 1 - /l 
-r 2 s > f 2 



V 



r s > 1 
r 2 s > 1 



h 
h 



v 



-r x s > fx 
r 2 s>l- /j 



(5.2) 



with s > 0, and the system (I2.8P of Theorem 1 2 . 1 1 b ecomes 



a 


+r l v\ 




+r 2 wi 


> o 


a 


-r 1 v 2 




+r 2 w 2 


> o 


a 


— r 1 ^ 




— r 2 W3 


> o 


a 






— r 2 W4 


> o 




+ f\Vl 




+f2Wl 


= 1 






)v 2 




= 1 






)v 3 


+ (1 " h)w 3 


= 1 








+(1 - h)wt 


= 1 


Vi, 


lOi > i 


e{i...4}. 





(5.3) 



where Vi,Wi, 



1, . . . , 4 stand for un,Ui2, i 



l,...,t = 2 q (since q = 2, t = 2 q = 4). 
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By Proposition 12.21 the cuts generated by the CGLP with constraint set (|5.3p and 
objective function mm pa for some p > have the form as > 1, where 

ctj = max{a], aj , aj , a j } 

with 

a] = —TjVl —TjlVl 



+r}v 2 -r 2 w 2 



I T V (5-4) 
- +rp 3 +r*w 3 



3 3 '3 

oq = -r l -Vi +t 2 -wa. 

As discussed in Section El a cut produced by the CGLP can be viewed as an 
intersection cut derived from a parametric cross-polytope or octahedron. For given 
v,w, we call the polyhedron 

Pocta,(v, w) = {(xi , x 2 ) G M 2 : -vi_xi - wix 2 < ; 

+v 2 x\ - w 2 x 2 < v 2 ; 
+w 3 xi + w 3 x 2 <v 3 +w 3 ; 

— V4X1 + U)4 1 X 2 < Wi } 

the (v,w) -parametric octahedron. 

If Vi = or Wi = for some i £ {1, . . . , 4} the z-th facet of P cta is parallel to one of 
the coordinate axes. If u,, vii > then the i-th facet of P cta is tilted (note that since 
we use the normalization f3 = 1, vt and Wi cannot both be 0). Varying the parameters 
v, w, the (v, ii?)-parametric octahedron produces different configurations according to 
the non-zero components of v,w. Depending on the values taken by the parameters, 
Poct&(v,w) may be a quadrilateral (i.e. a full-fledged octahedron in M 2 ), a triangle, 
or an infinite strip. In the rest of the section we refer to these configurations using 
the short reference indicated in parenthesis. It can easily be verified that the value- 
configurations of the parameters i?j , wi which give rise to maximal convex sets are the 
following: 

• (S) If exactly 4 components of (v, w) are positive, P cta is the vertical strip {x G 
R 2 : < x\ < 1} if V{ > 0, i = 1, . . . , 4; or the horizontal strip {x € K 2 : < x 2 < 
1} if Wi > 0,i = 1,...,4 (see figure |3 (a) j[3(b)| ). 



(Ta) If exactly 5 components of (v, w) are positive, P cta is a triangle with 1 tilted 
face (type A) (by "tilted" we mean a face that is not parallel to any of the two 



axes). Figure 3(c) illustrates the case with v\, w 2 , V3, 103, V4 > 0; wi, v 2 , W4 = 0. 
When in addition Vi = Wi for some i G {1, ... ,4} P cta becomes a triangle with 
vertices (0,0); (2,0); (0,2) or one of the other three configurations symmetric to 
this one. This corresponds to what is called a triangle of type 1 in [19J. In the 
general case Ta corresponds to a triangle of type 2 in |19j . 

(Tb) If exactly 6 components of (v, w) are positive, -P cta is a triangle with 2 
tilted faces (type B). Figure [3(d) | illustrates the case with vi, Wi, v 2 , w 2 , W3, W4 > 
0;f3,^4 = 0. This configuration corresponds to a triangle of type 2 in |19j . 

(Q) If all 8 components of (v,w) are positive, P cta is a quadrilateral. See Figure 



16 



The case with 7 components of (v, w) positive does not correspond to a maximal 
parametric octahedron, therefore we do not need to consider it. Suppose all the com- 
ponents are positive except for v\ which is 0. The facet of P cta corresponding to (0, 0) 
is horizontal and goes through the point (1,0). Is not hard to see that setting v 2 = 
we enlarge the set defined by the parametric octahedron. 




(a) 4 non-zeros - vertical strip (b) 4 non-zeros - horizontal strip (c) 5 non-zeros - triangle of type A 




(d) 6 non-zeros - triangle of type B (e) 8 non-zeros - quadrilateral 



Figure 3: Configurations of the parametric octahedron for the MIP case 
For a cut X^jeJ a j s j — 1 Andersen et al. [lj introduce the set 



L, 



x£M 2 :(x,s)eP L A^ ajSj < 1 



Clearly, L a C P oc t a (v, w), and the inclusion is often strict. 
Example. 

In PP, Andersen et al. considered the two rows instance 

\ +2si +ls 2 -3s 3 



+ls 5 



(5.5) 



| +lsi +ls 2 +2s 3 -ls 4 -2b s , (5.6) 



X 2 - 2 

xi, x 2 € Z, s > 



We present the complete description of the disjunctive hull for (|5.6p . In order to do 
so we generated the CGLP of (15. 6p using the normalization constraint (3 = 1 and we 
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considered all feasible bases. The CGLP produces 5 different facets. For each of these 
we show the configuration of the parametric octahedron that yields the corresponding 
cut in terms of the v, w variables: 

1. Cut (T B ): 2si + 2s 2 + 4s 3 + s 4 + ^s 5 > 1 
v 1 = 2; v 2 = f ; v 3 = 0; v A = 

w\ = 1; w 2 = |; u> 3 = 2; u> 4 = 2 

2. Cut (T B ): Isi + |s 2 + #s 3 + |s4 + U 5 > 1 

Wl = |; ^2 = 0; 1U3 = 0; u> 4 = ^ 

3. Cut (T A ): fai + 2s 2 + 4s 3 + «4 + |s 5 > 1 
vi = 2; u 2 = | ; ^3 = 0; v 4 = 

uii = 1; = 0; = 2; u) 4 = 2 

4. Cut (5): fsi + |s 2 + 12s 3 + |s 5 > 1 
vi = 4; v 2 = |; «3 = f ! ^4 = 4 

u;i =0; W2 = 0; W3 = 0; u> 4 = 

5. Cut (Tb): 2si + 2s 2 + f s 3 + f s 4 + ^s 5 > 1 

«i = ¥; «2 = I; w 3 = 0; -u 4 = 

= 7; ^2 = 7; ^3 = 2; u; 4 = 2 

Of the 5 facets of Pjj, 3 are facets for Pj: cuts 1, 2 and 4. Note that cut 4 is a split 
cut and can be derived using only the tableau row corresponding to the variable x 2 . 
Cut 3 and 5 are facets of Pjj by Theorem 12.81 

The condition given in Theorem 14.11 for an inequality ax > 1, facet defining for the 
disjunctive hull, to also define a facet of the integer hull specializes for the case q = 2 
to the following. For each of the four vertices p 1 of K, p l must lie on the line segment 
between two intersection points of rays r 3 with the boundary of P cta- As discussed in 
Section [3l the inequalities ax > 1 can be generated in a subspace of < 2 q = 4 variables, 
and then lifted into the full space by using the multipliers Wj), i = 1, . . . , 4. In [25] 
two algorithms were implemented for generating facets of the integer hull from Poctsu 
one for the case of a quadrilateral, the other for the case of triangles, both of them 
linear in \J\, the number of rays. 

Recently Dash et al. [18] have generalized the approach of [SI [6] , by considering 
more general 4-term disjunctions that give rise to what they call cross cuts and crooked 
cross cuts. They relate the closures of their cuts with the split closure and show, among 
others, that any 2 dimensional lattice free cut can be obtained as a crooked cross cut. 

6 Cut Strengthening 

Given a facet as > 1 of the disjunctive hull, if some non-basic variable Sj is required 
to be integral in the original problem formulation, then the cut can be strengthened. 
Let J\ be the index set of the integer-constrained variables Sj, and let J 2 = J \ J\. 

Lemma 6.1. If the disjunction 

-r l s>fi\ ( r l s>l-h \ fr l s>l-h\ f-r l s>h 
-r 2 s > h) V \-r 2 s >f 2 ) V V 2 s > 1 - h) V V r 2 s>l-f 2 
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where s > and Sj G Z, j G J\ C J, zs valid for Pj, then so is the disjunction obtained 
from \6. 1}) by replacing some or all r*-, i = 1,2, j G Ji, wt/i rj — m*-, for any m l j G Z, 
i = 1,2, j G J x . 



Proof. Suppose there exists i* G {1,2} and G J\ such that replacing r** with r** 
m** , where m** G Z, violates ()6. 1 j) . Then there exists a solution (x, s) G P/ with x G 2 
such that 

(_( r £ _ _ r)* Sj <f u ) A ((r£ - fh l 0s 3 , + ^ r}* Sj < 1 - ft 

holds. Rewriting this expression so as to bring together the terms in m** we get 



■2 



or 



-1 < m£ < 



contrary to the fact that both m** and Sj t are integer. □ 

Theorem 6.2. Given (v, w) > defining a parametric octahedron, the cut as > 1 can 
6e strengthened to as > 1 lozi/i coefficients aj,j G Ji given oy t/ie 3-variable mixed 
integer program 



mm a,- 



Q: ; 



—viuij — W\m 2 > —rjvi — r 2 w\ 
+V2m) - W2m 2 - > +r)v 2 - r 2 w 2 .„ 

-1-9 1- 9- (d.2) 

+t) 3 m^ + WzthTa > +r 1 f 3 + Hu> 3 

-1—9 1 — 9 - 

CKj —VilTlj + WiTUj > —TjV4 + rju^ 

iJipm 2 G Z. 

T/ie coefficients for j G J2 remain unchanged at aj = ctj as in Proposition \2.SX 

Proof. Validity of as > 1 follows from Lemma 16. 11 □ 

Theorem 6.3. The mixed integer program h6.2\) has an optimal solution (aj,m}j,m 2 ) 
satisfying G {|/<J> i = 1)2. 

Proof. Let (aj,m l j,m 2 ) be an optimal solution to the problem obtained from (|6.2p by 
adding the constraint m 1 - G l^j'l}- We wm snow that this solution cannot be 

improved by replacing rhpfh 2 with any other pair of integers. 

Consider the linear programming relaxation of (j6.2|) , which asks for minimizing the 
maximum of four linear functions. This is a piece-wise linear convex programming 

problem whose minimum is attained for m l - = rj, i = 1, 2, yielding aj = a j = = 

af = 0. From the convexity of the objective function a(mj,m 2 ) it follows that the 
integer optimum occurs at one of the points (m l -,m 2 ) G {(L r }J> L r |J)i (L r jJ) \ r j~\)> 
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([rj], L r f J)> (r r jl> r r |l)}- F° r suppose the optimum were to occur at some other point, 
say (mj,m|), where fh 1 - = [Y-] and m| = [r|] + dj for some dj > 0. Then 



aflr/l.r?) < a( [r/] , frj] ), 

hence 

a( [>/ 1 , frjl ) > Aa( \r}] , r?) + (1 - A)a( [r/] , [r|] + d,-) for < A < 1, 



i.e. the value of the minimum at a point which lies on the line between (\rj and 
[V|] + dj) is larger than a convex combination of the values of the minimum 
at the endpoints of the line, contrary to the assumption that a(m^,mj) is a convex 
function. □ 

The operation of replacing r*- by r*- — m*- for some m*- E Z, i = 1, 2, in the expression 
for a, is called the modularization of r 1 -, or more generally, the modularization of the 
cut ax > 1. Using m\ € {L r ]J> \ r j~\} is called the standard modularization. It can 
be shown (see below) that the mixed integer program (|6.2p attains its optimum for a 
standard modularization. 

Lemma 6.4. There exists a standard modularization f of the ray r such that 

0</i + f*<l, »€{1,2} (6.3) 
i.e. £/ie point (/ + f) belongs to K . 

Proof. If fi + r % — [r l \ < 1 then let m' = [?~*J • Note that the condition fi + r % — \ r % \ > 
follows since < fi < 1 and r* — |_r-* J > 0. Otherwise (/j + r l — \ r l \ > 1) let m l = \r l ~\ 
and from /j < 1 and r 4 — \r % \ < 1 we get < fi + r l — [r 1 ] — 1 = /, + r % — \r l ] < 1. 

For = 1, ... ,4, let a| be obtained from ag of (|5.4|) by substituting fj for r*-, 
i = 1,2. One can show that each is the convex combination of one of the expressions 

_ 1 _ 1 _ 2 — 2 

-j^ 2 - or j-f^- with one of the expressions or yzr^- To be specific, we have 
Lemma 6.5. 





-f 1 

= A ^ 




-r 2 

. 


with Ai 


= 






f 1 

= ^4 


+ (i 


-r 2 

-A 2 )^ , 


with A2 


= -v 2 (l 


-h) 


a) 


r 1 

= 


+ (i 


f 2 


wi/j A3 


= us(l 


-h) 


3 


-r 1 

= ^ A ~ft 




f 2 

A4) !_ /2 , 


wif/i A4 


= U4/1 





Proof. By substituting for the A&, k = 1, . . . , 4, we get the corresponding expressions 
for a*?. □ 

Theorem 6.6. T/ie strengthened cut as > 1 satisfies < ay < 1, j € J%. 
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Proof. Since Vk,Wk > for all k, we have a^ > for at least one of the four k, hence 
a.j > 0. Let (a,fh 1 ,fh 2 ) be an optimal solution to (|6.2p . Let f % = r % — fh % , i = 1,2, 
where m 1 € {L^J; I" 7 "*]} i = 1, 2. There are four cases: 
Case 1. m 1 = |_r*J , i = 1,2. Then f* = r* — m ! > 0, i = 1, 2, and 

a 1 = — f 1 ^! — f 2 tDi < 0. 

a 2 = f 1 ^ - r 2 w 2 < f x /(l - /i) (from (jO]l ). From ^/(l - fx) < 1, hence a 2 < 1. 

a 3 = f 1 ^ + r 2 wz = A3 f 1 /(l — fx) + (1 - As)f 2 /(1 - / 2 ), with A3 = £3(1 — /1) (from Lemma [675]) . 
But from Lemma l6.4[ f l /(l — /j) < 1, i = 1, 2, hence a 3 < 1. 

a 4 = — f V4 + f 2 u)4 < f 2 /(l — f%) < 1 (from [673]) . hence a 4 < 1. 

The remaining three cases, namely {m l ,m 2 ) = (|~f ], L^ 2 J)) {fh l ,fh 2 ) = (L^ 1 ], f^ 2 ])) 
and m ! = ([f ], [r 2 ]), i = 1,2, are similar. □ 

A way to further strengthen these cuts consists in the following three-step proce- 
dure: 

1. Apply standard modularization to each of the two rows from which the cut is 
generated (i.e. replace the ray r'- by rj — |_r*- J if r* > and by rj — [r*-] if r*- < 0, 
i = 1,2, jeJx). 

2. Generate a cut ax > 1 from the two modularized rows. 

3. Modularize the resulting cut to obtain the strengthened cut ax > 1. 

Yet another way to use the integrality of the variables Sj, j G Jx, is to apply the 
monoidal cut strengthening procedure of [9]. For cuts generated from a disjunction 
of the form f[6. 1|) . this procedure involves the use of lower bounds on the expressions 
on the lefthand side of each inequality. While these bounds are readily available and 
quite tight in the case when x±,X2 € {0, 1}, they can be weak in the general case of 
xx, X2 G Z. We therefore defer the dicussion of this procedure until the section on the 
0-1 disjunctive hull. 

7 The 0-1 Disjunctive Hull 

We now consider the 0-1 disjunctive hull Pp for q = 2, i.e. we work with Pqi = 
{(x,s) G {0, l} 2 x ]Rl J l : (x,s) e P L } where P L is given in ([571]). The CGLP that 
produces the facets of P^ is the linear program with the constraint set of Theorem [275] 
In addition to the four configurations of the parametric octahedron for the MIP CGLP 
given in Section [5j when v, w are unrestricted in sign some additional configurations 
are possible: (a) triangles with each face containing exactly one vertex of K, which we 
call triangles of type C (2b); and (b) cones, designated as (C). 

Note that our triangles of type C are similar to the class of triangles of type 3 for 
cuts for mixed integer programs described in [19] . The difference between these classes 
is that on the one hand, the three integer points contained in the faces of triangles of 
type 3 defined in [19] need not be vertices of K; on the other hand, our triangles of type 
C may also contain (non-0-1) integer points, positive or negative, in their interior. The 
presence among the parametric octahedra of unbounded ones, namely cones, implies 
that the cuts as > 1 of this class may have coefficients ay < 0. 
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As we did in Section [5l we give a classification of the parametric cross-polytopes 
that correspond to disjunctive hull facets for the 0-1 case (i.e. facets of Pp). Let k% € 
{1, . . . , 4} be the index of any vertex of K. We denote by k 2 , ks, k A the indices of the 
vertices of K that follow k\ in counter-clockwise order. The following configurations, 
in addition to those for facets of Pjj, are exhaustive when considering every value for 
k\ € {1,... ,4} (mod 4) and swapping Vi with Wi. In each case, the shape of P cta 
is determined by a strict subset of the four pairs (vi,Wi), the remaining pairs being 
inactive. 

• {Tex) v kl > 0,w kl < 0; v k2 ,w k2 > 0; v k , 3 > Q,w k3 > 0, (v k4 ,w k4 > 0). P octa is a 
triangle of type C with all its vertices outside the cube K. The face corresponding 



to A)4 is inactive. See Figure 4(a) 



{Tc2) v kl > 0,w kl < 0; v k , 2 ,w k , 2 > 0; v ks < 0,w ka > 0, {v ki ,w ki > 0). P oc ta is a 
triangle of type C with one vertex in the cube K. The face corresponding to k A 
is inactive. See Figure |4(b)| 

{Ca) v k2 ,v ks > 0; w k2 = w ks = 0, v ki > 0,w k4 < 0, (v kl ,w kl > 0). P octa is a cone 
with one face containing two adjacent vertices of K, the other face containing one 



vertex of K. The face corresponding to k\ is inactive. See Figure 4(c) 

{Cb) v kl < 0,w kl > 0; v k3 > 0,w k3 < 0, v k4 ,w k4 > 0, (v k2 ,w k2 > 0). P oc ta is 
a cone with one face containing two nonadjacent vertices of K, the other face 
containing one vertex of K. The face corresponding to k% is inactive. See Figure 



4(d) 



{Cc) v kl > 0,w kl < 0; v k2 ,w k2 > 0; (v k3 < 0,w k3 > 0), (v k4 ,w k4 > 0). P oc ta is 
a cone with each face containing one vertex of K. The faces corresponding to k% 



and &4 are inactive. See Figure 4(e) 



{Car) v kl > 0,w kl < 0; v k2 ,w k2 > 0; v k3 ,w ka > 0; (v k4 ,w k4 > 0). P oc ta is a 
truncated cone with each face containing one vertex of K. The face corresponding 
to &)4 is inactive. 

(5) v kl < 0,w kl > 0; v k3 > Q,w k3 < 0, (v k2 ,w k2 > 0,v klk ,w k4 > 0). P octa is a 
tilted strip, each side of which contains one vertex of K. The faces corresponding 



to the remaining two vertices are inactive. See Figure 4(g) 

• {St) v kl ,w kl > 0; v k2 ,w k2 > 0; v k3 > 0, w k . 3 < 0; {v k4 ,w k4 > 0). P octa is a 
truncated (tilted) strip, each side of which contains a vertex of K. The face 
corresponding to k^ is inactive. 

Example Consider the Andersen et al. pQ instance, amended with the condition 
Xi € {0, 1}, i G {1, 2}: 

xi =\ +2s\ +ls2 -3s 3 +ls 5 

x 2 = \ +lsi +ls 2 +2s 3 -ls A -2s 5 (7.1) 

xi,x 2 e {0,1}, s>o. 

In section Owe listed the 5 cuts defining the facets of the disjunctive hull for this 
example, without the 0-1 condition. Using the stronger disjunction expressing the 0-1 
condition we obtain the following 12 cuts that define the facets of convP^. 
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(a) triangle of type C with (b) triangle of type C with (c) cone with one face con- 
all vertices outside K one vertex inside K taining two adjacent ver- 
tices of K 




(d) cone with one face con- (e) cone with each face (f ) truncated cone with 
taining two nonadjacent containing one vertex of K each face containing one 
vertices of K vertex of K 




(g) tilted strip (h) truncated tilted strip 



Figure 4: Additional configurations of the parametric octahedron for the 0-1 case 
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1. Cut (type S): 2.667si + 1.333s 2 + 12s 3 + Os 4 + 1.333s 5 > 1 
v\ =4; v 2 = 1.333; v 3 = 1.333; v 4 = 4 

w\ = 0; w 2 = 0; u> 3 = 0; w 4 = 

2. Cut (type T B ): 2.667si + 1.333s 2 + 4.889s 3 + 0.8889s 4 + 1.333s 5 > 1 
V! = 2.222; v 2 = 1.333; v 3 = 1.333; v 4 = 0.4444 

Wl = 0.8889; w 2 = 0; w 3 = 0; w 4 = 1.778 

3. Cut (type T B ): 2si + 2s 2 + 4s 3 + ls 4 + 1.714s 5 > 1 
Wl = 2; v 2 = 1.143; v 3 = 0; t> 4 = 

Wl = 1; u> 2 = 0.2857; w 3 = 2; w 4 = 2 

4. Cut (type 2.947si + 1.053s 2 + 5.263s 3 + 0.8421s 4 + 3.579s 5 > 1 
vi = 2.316; v 2 = 0.7719; v 3 = 1.895; v A = 0.6316 

wi = 0.8421; w 2 = 0.8421; w 3 = -0.8421; w 4 = 1.684 

5. Cut (type T C i): 1.63«i + 2.37s 2 + 8.444s 3 + 0.4444s 4 + 1.926s 5 > 1 
«i = 3.111; v 2 = 1.037; u 3 = -0.7407; v 4 = 2.222 

Wl = 0.4444; w 2 = 0.4444; u> 3 = 3.111; w 4 = 0.8889 

6. Cut (type T C2 ): 4.364si + 2.545s 2 + 3.273s 3 + 1.091s 4 + 0.3636s 5 > 1 
vi = 1.818; v 2 = 1.818; v 3 = 1.818; v 4 = -0.3636 

wi = 1.091; w 2 = -0.7273; w 3 = 0.7273; w 4 = 2.182 

7. Cut (type T C2 ): 3.765si + 3.059s 2 + 2.588s 3 + 1.176s 4 + 0.7059s 5 > 1 
vi = 1.647; v 2 = 1.647; v 3 = 0.7059; v 4 = -0.7059 

wi = 1.176; w 2 = -0.4706; w 3 = 2.353; w 4 = 2.353 

8. Cut (type C A )- 12si + 8s 2 + 12s 3 + 0s 4 - 4s 5 > 1 
vi =4; u 2 = 4; f 3 = 4; v 4 = A 

wi = 0; it; 2 = —4; w 3 = 4; w 4 = 

9. Cut (type C B ): 32s x + 20s 2 - 20s 3 + 4s 4 + 12s 5 > 1 
v\ = —4; v 2 =4; v 3 = 4; v 4 = —12 

wi =4; it; 2 = 4; u> 3 = —4; w 4 = 8 

10. Cut (type C B ): 12s 1 + 8s 2 + 44s 3 - 4s 4 - 4s 5 > 1 
vi = 12; v 2 = 4; w 3 = 4; v 4 = -4 

1 = —4; w 2 = —4; w 3 = 4; w 4 = A 

11. Cut (type C c ): -2si + 6s 2 + 52s 3 + 2s 4 + 4s 5 > 1 
«i = 0; w 2 = 0; w 3 = —8; v 4 = 8 

Wl = 2; u/ 2 = 2; tf 3 = 14; w 4 = —2 

12. Cut (type C B )- 8«i - 4s 2 + 12s 3 + 16s 4 + 44s 5 > 1 
Wl =4; t> 2 = -9.333; t> 3 = 12; v 4 = 4 

«;i = 0; u> 2 = 16; w 3 = —16; u; 4 = 

The above list of 12 cuts includes 3 of the 5 cuts defining facets of convPo, namely 
1, 2 and 4, which appear on our list in position 3, 2 and 1, respectively. The remaining 

2 facets of convPo, given by cuts 3 and 5, are redundant for convP^; namely, cut 

3 is a convex combination of cuts 2, 3, 6 and 7 on our list, while cut 5 is a convex 
combination of cuts 1 and 5 on our list. 

The number of facets of convP^ substantially exceeds the number of facets of 
convPf). In order to assess the impact of the two sets of cuts, we computed the 
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average integrality gap for 1,000 randomly generated objective functions. Adding the 
5 cuts valid for the 2-row MIP reduces this gap by 77%; while adding the additional 
cuts valid for the 0-1 case reduces 100% of the gap. 

Next we discuss the strengthening of valid cuts for Pp when some variables Sj are 
integer-constrained. Let J\ be the index set of such variables. 

First of all, we observe that the standard modularization procedure described in 
Theorem 16.21 for strengthening cuts for Pjj is not valid in the case of cuts for Pp. 
Indeed, Lemma 16 . 1 1 which underlies the correctness of the procedure in the case of Pd, 
is no longer valid in the case of P^: if the disjunction (|6.ip is modified by replacing 
every inequality with equality, then it is no longer equivalent to the disjunction obtained 
by replacing r l - with r l - — m 1 -. Instead, we will use a different modularization, known 
in the literature under the name of monoidal strengthening [9]. 

Consider a disjunction of the form 



\J(A k 



x>og), A k = (a*), j G J, 4eR m ,jGJU{0}, (7.2) 



keQ 

and the valid cut ax > 1, where 



max{0 fc aj/0 fc ag} (7.3) 



for some 9 k G M+ , k G Q. 

Suppose now that for each A k x, k G Q, a lower bound 6q < is known, i.e. 
A k x >b^,ke Q. 

Theorem 7.1. Let M := {m G : J2keQ mk ^ °}- V x 3 G j G Ji, then the cut 
ax > 1 can be strengthened to ax > 1, where 

aj = nun max { (e k a k + m k 9 k {a k - /6 k a k ) } j G J x (7.4) 
and ctj = aj for j G J\J\. 

Proof. See [9j or @]. □ 



We will now apply this Theorem to our case, first with the disjunction (|6.ip . then 
with the stronger disjunction defining Pp. Let ax > 1 be an inequality implied by 
the disjunction (|6.ip . i.e. a valid inequality for Pp, and let's assume that x% G {0, 1}, 
i = 1,2. It is not hard to see that a lower bound on the lefthand side of each of 
the 8 inequalities that occur in (|6.ip is obtained by subtracting 1 from the righthand 
side. This means that if <2q denotes the righthand side and 6g the lower bound on the 
lefthand side of the k-th term, then Oq — 6g = Q) . 

Now the cut from the disjunction (|6.ip is as > 1, where 



and 



ay = 


max {a k \, 
fce{l,...,4} 




-fjvi - r 2 -w\ 




rjv 2 - r|u) 2 


a| = 


rjV S + rjws 




-TjV^ + r|w} 4 



(7.5) 
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To apply the theorem to this case, notice that 9 k = (vk,Wk) and 6 k a k = a k , 9 = 1 
for k = 1, ... ,4. 

Corollary 7.2. Let x u x 2 G {0, 1}, and I={m£Z 4 : Yl=i mk > °}- T/ien « s > 1 
is a va/ie? cut for Pd, with aj = max/ cg { 1) ... ir }{a k }, and 



n k 
j 



min max {a k + m k (vk + w^)} j G Ji 
a^eMfce{i,...,4i 



m^GM fc6{l,...,4} 



Proof. Denoting 
i /7i 



fl-1 

/2-1 



; "0 



1 - / M a 3 



-/i 



! "0 



1-/2 

-/l 
"/2 



J 2 - K 

it is easy to see that for k = 1, . . . , 4, 

(ujfc, u) fc )(ao - &o) = ^fc + ^fc- 



j € J \ Ji 



, "0 



/l 



6 » = ,7 -A 1 )' 



(7.6) 



□ 



We now turn to strengthening a valid inequality for P^, the set defined by the 
disjunction (|6.1F ). obtained from (|6.ip by replacing each inequality with equality. In 
this case the cut from (16. IT ) is ax > 1, where Sy = max^g^ 4 } 5^ and the 5^ are given 
by the same expressions ()7.5f) as , with the important difference that the parameters 
(vk,Wk) are unrestricted in sign. However, from the normalization constraints (|5.3p it 
follows that for any k G {1, . . . , 4}, at most one member of the pair {vk^Wk) can be 
negative. 

In order to derive the lower bounds 6g required by Theorem 17. 1\ the best way is to 
represent each equation of (j6.1T ) as a pair of inequalities; i.e. the first term of (j6.1T ) 
is restated as 

/ _ r i s > h \ 



V 



2 

-r s 



2 

r s 



> 



> 



> 



-fx 
h 

-h ) 



(7.7) 



and so on. Denoting the corresponding parameters or multipliers by v' k ,v' k \w' k ,w' k ' for 
k = 1, ... ,4, we see that since at most one of the pairs of inequalities corresponding 
to an equation can be active in any given solution, at most one member of each pair 
(v' k ,v k ) can be positive, and the same holds for each pair (w' k ,w k ). Furthermore, it 
becomes clear that if in the equality formulation (j6.1T ) a parameter, say v\, takes on 
a negative value v\ < in a solution, this corresponds to the fact that the member of 
the pair of inequalities corresponding to the equation associated with v\ that is active, 
is the one with <, i.e. with the inequality reversed. 
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Corollary 7.3. Let M = {m E Z 4 : Ylt=i mk ^ °l- T/ien 

as > 1 is a va/zd cut for 

Pp, with ctj = maxfcgj! )4 }{a^} ; and 



o 




max {a* + m^(>£ + io^)}, j € J\ 



j€J\Jl 



where v,' 



maxjufc, 0} and wt = maxj^ , 0} . 



Proof. If, using the inequality formulation (|7.7p of the disjunction (|6.1 
the righthand sides of the four terms by 



«o 



/ 


fl 


\ 




fl" 


-h 


\ 




( 1- 


\ 




-fl 






h 


- 1 






/i 


- 1 




h 




, al = 




/ 2 






1 - 


-/ 2 


V 


-h 


) 




\ ~ 


h 


) 




u 


-1 / 



~ 1 
"o 



we denote 



/l \ 

~f2 
1-/2 

V /a - 1 y 



(7-8) 

then the lower bounds on the expressions on the lefthand sides of the inequalities are 
no longer equal to the righthand side minus 1. Instead, we have the following situation: 



/ fl ~ 1 \ 

-fi 

/2-1 

V -h ) 



As a consequence, 



/i-l 

/2-1 

V -h J 



( ~h \ 
/i-i 
-h 
V h - 1 J 



( h - 1 \ 
-h 
-h 
V h - 1 J 



(7.9) 



n K — h K 



f 1 \ 


1 

V o y 



for k = 1, 



Thus, if we denote := max{t>fc,0}, = max{w?fc,0}, we have (v£, w~£)(o,q, 6q) 



(vjf,u£), fc 



, 4, and the expression for a*? follows. 



□ 



Finding the optimal rrij G M requires a small (single digit) number of comparisons. 
While [9] and [I] give simple procedures for the case of a general disjunction, the 
optimal rrij of Corollary 17.31 for a given j 6 Jx can be found as follows: 

• Start with m k - = for all k and apply the 
Iterative Step. 

min/j a k and let m° 



Find af ax 



sponding values of m 
Set m™ ax = m 



max fc a^, af m 

k 

r 



max m mm corre _ 



m - 



+ t, where i is the smallest positive 



integer for which the identity of a™ ax changes. 
If the value of max^ a* has not been reduced, stop with m™ 



-- mfx—t+l, 

m r ? m + t — l, and mh unchanged for k ^ max, min. Otherwise repeat. 
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In the case where P cta is a triangle with each face containing exactly one vertex of 
K, the term of the disjunction (|6.1F ) corresponding to the vertex of K left outside the 
triangle plays no role in defining the cut, hence it can be dropped and the strengthening 
becomes simpler. This is even more true of the case of a cone, where only two terms 
of the disjunction are active. A particularly simple case is that of a "fixed shape" 
cone with apex at a vertex of K, and one face containing a side of K, the other 
face containing the diagonal of K. There are eight such cones, and every fractional 
(/lj/b) 7^ (55 2) (i- e - n °t ly in g on the diagonal of K) is strictly contained in four of 
them (see figure [6] in the next section). 

We will illustrate the monoidal strengthening procedure on the conic cuts obtainable 
from these disjunctions. Here is a couple of them: 



Each disjunction violated by the point (/1, $2) has positive righthand sides and gives rise 
to a valid cut as > 1, with coefficients ctj shown below, obtained by using multipliers 
normalized to yield a cut with a righthand side of 1: 



To apply the strengthening procedure, we note that for each of the 16 terms of 
the above 8 disjunctions, the lower bound on the lefthand side of the inequality is just 
1 unit less than the righthand side, hence the difference between the latter and the 
former is exactly 1. Further, the weights (vk,Wk) are normalized so that Vk + Wk = I, 
k = 1,2. The resulting strengthened coefficients for the above illustration are 



In this section we present computational experiments with cuts derived from fixed 
configurations of the parametric octahedron. We assess the strength of the cuts by 
analyzing the gap closed on instances from MIPLIB3_C_V2 |23j when used in com- 
bination with standard Gomory cuts. MIPLIB3_C_V2 is a collection of 68 instances 
by Margot which are slight variations of the standard MIPLIB3 [Mj and for which 
the validity of a candidate solution can be checked in finite precision arithmetic. We 
restricted the collection to a subset of 41 instances. The considered instances are such 



1. (-x 2 > 0) V (-xi + x 2 > 0) 

2. (x 2 > 1) V(-xi -x 2 > -1) 



or, after substituting /j + r l s for Xi, i = 1,2, 

1. {-r 2 s>f 2 )V{(-r l +r 2 ) S >h-f 2 ) 

2. {r 2 s > 1 - / 2 ) V ((-r 1 - r 2 )s > h + f 2 - 1). 







8 Computational Experiments 
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that they contain at least 2 binary variables fractional in the optimal LP solution and 
the cut generation procedure on each round takes less than 3600 seconds. 
We generated the following two families of cuts 

• Cuts from 4 Triangles Ta (shown in Figure [5]) whose vertices, expressed in terms 
of their x\,X2 coordinates, are: 

— (0,0); (2,0); (0,2) 

— (-1,0); (1,0); (1,2) 

— (0,-1); (2,1); (0,1) 

- (1,-1); (1,1); (-1,1) 

• Cuts from 4 of the 8 cones of type Ca (shown in Figure [6]): 

— apex at (0,0) and rays (1,0), (1,1) 

— apex at (0, 0) and rays (0, 1), (1, 1) 

— apex at (0,1) and rays (1,0), (1,-1) 

— apex at (0, 1) and rays (0, —1), (1, —1) 

— apex at (1, 1) and rays (—1, 0), (—1, —1) 

— apex at (1, 1) and rays (0, —1), (—1, —1) 

— apex at (1,0) and rays (— 1, 0), (— 1, 1) 

— apex at (1,0) and rays (0,1), (—1,1) 

The reason we only used 4 of these 8 cones is that every (/i, /2)-pair is contained in 4 
of these 8 cones. 




-10 12 3 -2-1012 




-10 12 3 -2-1012 



Figure 5: Fixed shape Triangles Ta 

For each instance, we first solved the linear programming relaxation and generated 
a round of Gomory mixed integer (GMI) cuts, a round being one cut from every row of 
the optimal simplex tableau associated with a binary basic variable with a fractional 
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1 



2 



3 



4 




Figure 6: Fixed shape cones Ca 



value. We then generated from each pair of rows with at least one fractional binary 
basic variable either (a) all cuts from the 4 triangles Ta, or (b) all cuts from 4 of the 
8 cones Ca or both, and strengthened them via standard modularization (in case (a)) 
or monoidal strengthening (in case (b)). 

We call this cut generating cycle a round. At the end of each round, we reoptimized 
the resulting linear program and removed all cuts that were not tight at the optimum. 
We generated up to 5 rounds of cuts for each instance. A statement of our routine 
follows. 
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Cut Generating Procedure (r, /) 

1 Solve LP relaxation P 

2 for k ^— 1 up to 5 



3 do 

4 Initialize cut collection C ^— empty 

5 for each binary basic Xi fractional in the 
current solution 

6 do 

7 Compute Gomory cut Gi 

8 C ^CUd 

9 for each binary basic pair Xi , Xj with at 
least one fractional in the current solution 

10 do 

11 if GenerateTriangles==true 

12 then generate the cuts T^, . . . , from each of the 4 triangles of 

Ta that contain the fractional solution in their interior 

13 if StrengthenCuts==true 

14 then strengthen the cuts Th, . . . , via standard modularization 

15 c< cut* v;- 

16 if GenerateCones==true 

17 then generate the cuts Kfj, . . . , Kf- from each of the 8 cones Ca that 

contain the fractional solution in their interior 

18 if StrengthenCuts==true 

19 then strengthen the cuts K^, . . . , Kfj via monoidal cut 

strengthening 

20 C i-CUK!,...,^ 

_, • if. 

21 Resolve P and get new solution x with value opt 

22 Remove from P the cuts in C that are not tight at x k 



The Gomory mixed integer (GMI) cut generator we used is the CglGomory routine 
of the Cgl package of COIN-OR |14j . Tables [87018,41 summarize the results of our 
experiments with these cuts. Table 18.11 shows the outcome of applying all three types 
of cuts in the above described manner, with strengthening, for one round. Column 
1 lists the 41 test instances mentioned above. Column 2 shows the percentage of the 
integrality gap closed by one round of GMI cuts, while the next two columns show the 
number of cuts generated and added to the LP relaxation, along with the number of 
cuts deleted after reoptimization as nonbinding. The next three columns show the same 
data (i.e. percentage of gap closed and number of cuts added, respectively deleted) 
after generating a cut from each of the 4 triangles Ta associated with every pair of 
basic 0-1 variables with at least one fractional member, and a cut from each of the 8 
cones associated with every such pair, provided the cone contains such a pair in its 
interior. Finally, the last column shows the percentage improvement in the integrality 
gap closed by all three types of cuts versus the GMI cuts alone. 

As the table shows, the integrality gap closed, which is 19.49% in the case of the 
GMI cuts, reaches 29.06% when the two remaining types of cuts are added, an increase 
of 49.14%. The number of triangle cuts and conical cuts generated is of course much 
larger than that of GMI cuts. While the latter is bounded by the number of basic 
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0-1 variables fractional at the optimum, in case of the other two types of cuts this 
number gets multiplied by 8 times the number of basic 0-1 variables, fractional or not. 
From the table it is clear that after reoptimization few of the added cuts remain active 
(about 2%), while the rest get removed. The table also reveals marked differences in 
the impact of the 2-row cuts on different instances, from impact in about 40% of the 
instances, to a more than 7-fold increase of the gap closed in the highest-impact case. 

Tables I8.2H8.3I show the effect of using only triangle cuts or only conic cuts on top 
of the GMI cuts. Clearly, the joint effect of using both types of cuts is substantially 
stronger than is the case with a single type. 

Finally, Table 18.41 shows the effect of generating both types of 2-row cuts on top of 
GMI cuts, as in Table [8TTT but this time for 5 rounds instead of just 1. The improvement 
in gap closing keeps growing after every round. At the end of the 5 rounds, the gap 
closed is 38.78%, roughly twice as large as the 19.49% gap closure obtained by 1 round 
of GMI cuts. 
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Table 8.1: GMI cuts + Triangle cuts + Conic cuts, all strengthened, 1 round 





GMI 


GMI+ 


T A + C A , 


strengthened 






Gap 


Cuts 


Cuts 


Gap 


Cuts 


Cuts 






closed 


added 


deleted 


closed 


added 


deleted 


Improvement 


Instance 


% 


# 


# 


% 


# 


# 


% 


air03 


100 


36 


10 


100 


127 '44 


12667 


0.00 


cap6000 


41.65 


2 


1 


41.65 


1360 


1359 


0.00 


danoint 


0.26 


24 


13 


0.26 


24 


13 


0.00 


dcmulti 


45.75 


49 


12 


45.86 


392 


362 


0.24 


egout 


21.84 


16 





60.91 


3461 


3423 


178.89 


enigma 


100 


6 


5 


100 


341 


340 


0.00 


fiber 


53.32 


39 


24 


64.97 


22944 


22881 


21.85 


fixnct3 


6.62 


6 





56.19 


2925 


2827 


748.79 


fixnet4 


4.79 


6 





13.02 


2829 


2731 


171.82 


fixnet6 


3.98 


6 





13.03 


2502 


2388 


227.39 


khb05250 


74.91 


19 





84.21 


1435 


1408 


12.41 


11521av 











26.68 


16774 


16684 


0.00 


lseu 


55.94 


12 


7 


56.59 


567 


562 


1.16 


marksharel 





6 


3 





124 


110 


0.00 


marksharc2 





7 


3 





173 


147 


0.00 


mas74 


6.52 


9 





7.57 


511 


485 


16.10 


mas76 


6.36 


9 


1 


7.7 


433 


412 


21.07 


misc03 


8.62 


20 


17 


8.62 


2239 


2231 


0.00 


misc06 


26.17 


8 





26.17 


8 





0.00 


misc07 





28 


25 


0.72 


4177 


4171 


0.00 


mod008 


20.1 


4 


1 


20.27 


96 


90 


0.85 


modOlO 











99.26 


14211 


14070 


0.00 


modOll 


11.44 


8 


1 


32.81 


855 


171 


186.80 


modglob 


13.32 


16 


2 


16.26 


137 


40 


22.07 


p0033 


12.6 


5 


1 


57.04 


174 


168 


352.70 


p0201 


16.89 


20 


13 


19.31 


1933 


1929 


14.33 


p0282 


3.47 


24 


17 


6.2 


2837 


2829 


78.67 


p0548 


3.06 


19 


2 


18.53 


5584 


5521 


505.56 


p2756 


0.21 


7 


1 


0.56 


908 


884 


166.67 


pkl 





15 


5 





22 


12 


0.00 


pp08a 




ou 


n 
u 


65.29 


7825 


7732 


20.24 




32.83 


40 





41 1 8 




9^1 3 


9^ A% 

AO .4:0 


qiu 


0.33 


36 


24 


0.33 


36 


24 


0.00 


rentacar 





2 








7 


5 


0.00 


rgn 


3.15 


17 


9 


3.15 


1341 


1331 


0.00 


setlch 


30.36 


125 


1 


44.51 


28347 


28111 


46.61 


stein27 





21 


17 





1257 


1254 


0.00 


stein45 





35 


28 





3584 


3575 


0.00 


swath 


8.18 


10 





11.79 


6917 


6895 


44.13 


vpml 


20.73 


12 





23.94 


1073 


1047 


15.48 


vpm2 


11.25 


27 


6 


16.96 


4345 


4296 


50.76 


Average 


19.49 


19.54 


6.07 


29.06 


3903.2 


3846.29 


49.14 
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Table 8.2: GMI cuts + Triangle cuts, strengthened, 1 round 





GMI 


GMI - 


- Ta, strengthened 






Gap 


Cuts 


Cuts 


Gap 


Cuts 


Cuts 






closed 


added 


deleted 


closed 


added 


deleted 


Improvement 


Instance 


% 


# 


# 


% 


# 


# 


% 


air03 


100 


36 


10 


100 


5646 


5569 


0.00 


cap6000 


41.65 


2 


1 


41.65 


364 


363 


0.00 


danoint 


0.26 


24 


13 


0.26 


24 


13 


0.00 


dcmulti 


45.75 


49 


12 


45.75 


119 


86 


0.00 


egout 


21.84 


16 





60.9 


2026 


1990 


178.85 


enigma 


100 


6 


5 


100 


143 


142 


0.00 


fiber 


53.32 


39 


24 


59.77 


9927 


9877 


12.10 


fixnet3 


6.62 


6 





47.01 


1571 


1479 


610.12 


fixnet4 


4.79 


6 





12.49 


1523 


1426 


160.75 


fixnet6 


3.98 


6 





12.03 


1352 


1240 


202.26 


khb05250 


74.91 


19 





84.21 


722 


698 


12.41 


11521av 











13.42 


7568 


7483 


0.00 


lseu 


55.94 


12 


7 


55.94 


291 


286 


0.00 


marksharel 





6 


3 





66 


56 


0.00 


marksharc2 





7 


3 





91 


69 


0.00 


mas74 


6.52 


9 





7.44 


261 


236 


14.11 


mas76 


6.36 


9 


1 


7.12 


225 


203 


11.95 


misc03 


8.62 


20 


17 


8.62 


1031 


1028 


0.00 


misc06 


26.17 


8 





26.17 


8 





0.00 


misc07 





28 


25 





1933 


1930 


0.00 


mod008 


20.1 


4 


1 


20.11 


48 


43 


0.05 


modOlO 











93.23 


5873 


5735 


0.00 


modOll 


11.44 


8 


1 


32.53 


464 


103 


184.35 


modglob 


13.32 


16 


2 


15.75 


107 


12 


18.24 


p0033 


12.6 


5 


1 


57.04 


85 


80 


352.70 


p0201 


16.89 


20 


13 


19.31 


1235 


1230 


14.33 


p0282 


3.47 


24 


17 


5.38 


1416 


1409 


55.04 


p0548 


3.06 


19 


2 


17.37 


2958 


2917 


467.65 


p2756 


0.21 


7 


1 


0.56 


546 


523 


166.67 


pkl 





15 


5 





18 


8 


0.00 


pp08a 




ou 


u 


64.89 


OOOZ 




19.50 




32.83 


40 





4(1 Q 


1132 


1058 


Zi4:.(JO 


qiu 


0.33 


36 


24 


0.33 


36 


24 


0.00 


rentacar 





2 








4 


2 


0.00 


rgn 


3.15 


17 


9 


3.15 


697 


689 


0.00 


setlch 


30.36 


125 


1 


44.3 


12412 


12154 


45.92 


stein27 





21 


17 





846 


843 


0.00 


stein45 





35 


28 





2395 


2386 


0.00 


swath 


8.18 


10 





11.79 


3358 


3338 


44.13 


vpml 


20.73 


12 





21.94 


480 


457 


5.84 


vpm2 


11.25 


27 


6 


16.02 


2133 


2094 


42.40 


Average 


19.49 


19.54 


6.07 


27.98 


1829.17 


1781.78 


43.61 
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Table 8.3: GMI cuts + Conic cuts, strengthened, 1 round 





GMI 


GMI H 


- Ca, strengthened 






Gap 


Cuts 


Cuts 


Gap 


Cuts 


Cuts 






closed 


added 


deleted 


closed 


added 


deleted 


Improvement 


Instance 


% 


# 


# 


% 


# 


# 


% 


air03 


100 


36 


10 


100 


7134 


7057 


0.00 


cap6000 


41.65 


2 


1 


41.65 


998 


997 


0.00 


danoint 


0.26 


24 


13 


0.26 


24 


13 


0.00 


dcmulti 


45.75 


49 


12 


45.86 


322 


292 


0.24 


egout 


21.84 


16 





35.15 


1451 


1422 


60.94 


enigma 


100 


6 


5 


100 


204 


203 


0.00 


fiber 


53.32 


39 


24 


62.94 


13056 


12970 


18.04 


fixnet3 


6.62 


6 





49.3 


1360 


1257 


644.71 


fixnet4 


4.79 


6 





10.41 


1312 


1210 


117.33 


fixnet6 


3.98 


6 





11.5 


1156 


1050 


188.94 


khb05250 


74.91 


19 





77.59 


732 


704 


3.58 


11521av 











26.68 


9206 


9118 


0.00 


lseu 


55.94 


12 


7 


56.59 


288 


283 


1.16 


marksharel 





6 


3 





64 


55 


0.00 


marksharc2 





7 


3 





89 


73 


0.00 


mas74 


6.52 


9 





7.15 


259 


237 


9.66 


mas76 


6.36 


9 


1 


7.62 


217 


190 


19.81 


misc03 


8.62 


20 


17 


8.62 


1228 


1220 


0.00 


misc06 


26.17 


8 





26.17 


8 





0.00 


misc07 





28 


25 


0.72 


2272 


2266 


0.00 


mod008 


20.1 


4 


1 


20.27 


52 


46 


0.85 


modOlO 











97.91 


8338 


8198 


0.00 


modOll 


11.44 


8 


1 


27.3 


399 


99 


138.64 


modglob 


13.32 


16 


2 


13.94 


46 


2 


4.65 


p0033 


12.6 


5 


1 


24.59 


94 


89 


95.16 


p0201 


16.89 


20 


13 


16.89 


718 


711 


0.00 


p0282 


3.47 


24 


17 


5.4 


1445 


1436 


55.62 


p0548 


3.06 


19 


2 


6.53 


2645 


2599 


113.40 


p2756 


0.21 


7 


1 


0.21 


369 


358 


0.00 


pkl 





15 


5 





19 


9 


0.00 


pp08a 


^A "} 


ou 


n 
u 


57.01 


401 3 


3Q97 


4.99 




32.83 


40 







1487 


1409 


4 9fi 


qiu 


0.33 


36 


24 


0.33 


36 


24 


0.00 


rentacar 





2 








5 


3 


0.00 


rgn 


3.15 


17 


9 


3.15 


661 


650 


0.00 


setlch 


30.36 


125 


1 


37.43 


16060 


15780 


23.29 


stein27 





21 


17 





432 


428 


0.00 


stein45 





35 


28 





1224 


1214 


0.00 


swath 


8.18 


10 





9.89 


3569 


3537 


20.90 


vpml 


20.73 


12 





22.73 


605 


575 


9.65 


vpm2 


11.25 


27 


6 


15.25 


2239 


2199 


35.56 


Average 


19.49 


19.54 


6.07 


25.88 


2093.56 


2046.59 


32.83 
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Table 8.4: GMI cuts + Triangle cuts + Conic cuts, all strengthened, 5 rounds 





GMI 


GMI + T A + C A , strengthened 






Gap 


Gap 






Closed 


Closed 


Improvement 


round 


% 


% 


% 


1 


19.49 


29.06 


49.14 


2 


24.94 


34.01 


36.38 


3 


27.87 


36.70 


31.67 


4 


29.57 


37.95 


28.31 


5 


30.48 


38.78 


27.24 
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